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Abstract 

We present a Monte Carlo generator for exclusive vector meson production in charged 
lepton-proton interactions, I + p — > I + p + V, based on a QCD leading logarithm model 
calculation according to which the cross section for this process is proportional to the 
square of the gluon momentum density in the proton. The generator can be used for 
both fixed target and collider kinematics. Each event is assigned a weight equal to its 
cross section, so that the independent variables can be generated according to convenient 
distributions, not necessarily identical to the physics ones. 



1 Introduction 

In this paper we present a Monte Carlo generator for exclusive production of vector mesons in 
the reaction 

e(n)+p^e(n)+p + V, (1) 



where an electron (e) or a muon (//) interacts with a proton (p), producing, in the final state, 
a vector meson V (p°, u>, 0, J/if)...). This reaction is also referred to as elastic vector meson 
production. 

Figure [I] shows the diagram of the process; tables [I] and |2| define the relevant kinematic 
variables. The incoming lepton radiates a virtual photon of mass — Q 2 that collides with the 
target proton, turning into a vector meson; the proton emerges intact from the interaction. 
Since the vector meson has the same quantum numbers as the photon, the object exchanged 
between the photon and the proton, a pomeron, carries the quantum numbers of the vacuum. 
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Figure 1: Elastic vector meson production. 
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four-momentum of the incident lepton 


k' 




four-momentum of the scattered lepton 


p 




four- momentum of the incident/target proton 


p' 




four-momentum of the scattered proton 


V 




four-momentum of the produced vector meson 


my 




mass of the vector meson 


m, p 




mass of the target proton 


q = 


k-k' 


four-momentum of the virtual photon 


-Q 


2 = q 2 


invariant mass squared of the virtual photon 


y = 


p-q/p-k 


fraction of the beam energy carried by the photon 






(in the proton rest frame) 


V = 


p ■ q/mp 


photon energy (in the proton rest frame) 


w -- 


= V(q + p) 2 


photon-proton centre of mass energy 




= V(k+P) 2 


electron-proton centre of mass energy 



Table 1: Definition of the relevant kinematic variables, part 1. 



t = 


--{p-p'f 


four-momentum squared exchanged at the proton vertex 


Pt 




transverse momentum of the vector meson with respect to 
the virtual photon direction 


z - 


= p-v/p-q 


fraction of the photon energy carried by the vector meson 
(in the proton rest frame) 


Xl 


= |p'l/|p| 


ratio between outgoing and incoming proton three-momenta 
in the laboratory frame 


p 

Pt 




tranverse momentum of the proton in laboratory frame 


X - 


Q 2 +m 2 v +p 2 t 
W 2 


fraction of the proton's momentum carried by the two-gluon 
system in the model of ref. |2l| 


f 


Q 2 +ra 2 v +p 2 t 


scale at which the gluon density is probed in the model of ref. 
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polar angle of the positive decay particle (e.g. tt + in the p° case) 
with respect to the outgoing p direction in the V rest frame 
(for two-body decays) 






angle between the decay plane and the meson production 
plane, which contains the virtual photon and the meson, in the 
V rest frame (for two-body decays) 


$ 




angle between the lepton scattering plane and the vector 
meson production plane (in the proton rest frame) 




= (p - $ 





Table 2: Definition of the relevant kinematic variables, part 2. 
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Pomerons were originally introduced to describe diffractive interactions between hadrons. 
Elastic vector meson production at small Q 2 can indeed be thought of as a hadronic process 
in the Vector Meson Dominance (VMD) [|TJ framework, in which the photon is assumed to 
fluctuate into a vector meson before interacting with the proton. Reaction (U) may thus allow 
the investigation of the properties of the pomeron, as well as of the hadronic features of the 
photon. 

Elastic vector meson production has been extensively studied in fixed target experiments in 
the photoproduction regime (Q 2 = 0) 0-PJ and at larger values of Q 2 10-0, for photon-proton 
centre of mass energies W < 20 GeV. With the advent of the electron-proton collider HERA 
reaction ([!]) has also been investigated at values of W ~ 100 GeV, both for Q 2 « [ 10 1 - [ 15 
and for non-zero virtualities f|TB| 



From the theoretical point of view, several models [0|-[23 offer a description of elastic 
vector meson production; all assume pomeron exchange, with the pomeron described as a gluon 
pair. Individual models differ in the way the gluons are treated, ranging from non-perturbative 
approaches to perturbative ones. 

The generator we wrote, hereafter referred to as DIPSI, is based on the model of ref. WW\ 
This model assumes that the exchanged virtual photon fluctuates into a qq pair which then 
interacts with a gluon ladder emitted by the incident proton (fig. |2]). The parameters of the 
model are the strong coupling constant a s , the two-gluon form factor of the proton and the 
gluon momentum density in the proton. Once these are chosen, the dependences on W and t, 
where t is the four-momentum squared transferred at the proton vertex, are fixed. 

DIPSI has been extensively used within the ZEUS experiment - operating at the electron- 
proton collider HERA at DESY - in particular for the determination of the detector acceptance 
for p° JlO], PH ) cf) [|lt| and J/ip photoproduction as well as for p° Jl6| and []Tj| production 



at large Q . 

The paper is organised as follows. In section ^| we give a short description of the model and 
in section ^ we present the program itself. The appendices contain a list of the files provided, a 
description of the input control cards as well as a list of the routines and functions used in the 
program. 



2 The model 



The model |2T|] is based on a perturbative QCD (pQCD) calculation of the 7 — > V process and 
assumes a non-relativistic wavefunction for the vector meson V. In the Born approximation, at 
small values of the parton fractional momentum x, the amplitude for the reaction 7^ — > Vp is 
given by the sum of those of the graphs shown in fig. ^| (without the "horizontal" gluons which 
form the rungs of the ladder), where the two t-channel gluons play the role of the pomeron (in 
the Born Low-Nussinov approximation ||25|| ). It is convenient to write the contribution of the 
upper quark loop in terms of the electronic width Tj e = T of the vector meson. In this way a 
large part of the 0(a s ) corrections is absorbed in the experimental value of T; here a s is the 
strong coupling constant. To take into account logarithmic corrections which contain (a s lnq 2 ) n 

1 The results of ref. lES lead to formulae with the same structure. 
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or (a s \nl/x) m terms (here q 2 denotes the scale relevant for the process), one has to sum up a 
set of more complicated Feynman diagrams. In the axial gauge, these diagrams have the ladder 
form (i.e. include all gluons shown in fig. |2|). At small t the lower part (that below the fermion 
loop) of the graphs of fig. ^| may be considered as the gluon momentum density in the proton 
xg(x,q 2 ). Indeed, for t = this lower part of the graph contains exactly the same Feynman 
diagrams (calculated for the same kinematics) as those involved in the definition of xg(x, q 2 ) in 
the Leading Logarithm Approximation (LLA) QCD treatment of deep inelastic scattering. 

The arguments x and q 2 of xg(x,q 2 ) take the values x = (Q 2 + rriy + p 2 )/W 2 and q 2 = 
(Q 2 + rriy + p 2 )/4, respectively. The value of x is larger than that of the Bjorken variable 
XBj — Q 2 /W 2 , as one has to transfer additional longitudinal momentum along the pomeron 
(i.e. the t-channel gluons) to put the vector meson on mass shell. 



The typical virtuality q was chosen equal to \ y Q 2 + m 2 / + p 2 as the t-channel gluons interact 
with one quark only, which, in the non-relativistic approximation used in the model, carries one 
half of the external momenta (i.e. those of the photon and of vector meson). 



At large energy the value of x is small and for zero angle production, when t 



^min — 



■x 2 m 



N 



gives: 



0, we can write the cross section in terms of T and g(x, q 2 ). An explicit calculation 



da T { , yp — > Vp) 
It 



a 2 {q 2 )Tm v 



3a 



xg(x,q 



7T X 



f(q 2 ,p 2 t' 



-i 2 



2q 2 (2q 2 -p 2 )\n(8q 2 /p 2 



F 2 N G (t) 



(2) 



where a = 1/137 is the electromagnetic coupling constant and the function f(q 2 ,p 2 ) is defined 
as follows: 



f(f,p 2 t)=^ 
f(f,Pt)=^ 



4q 2 - p 2 + p 2 



Pt + Po 



7-A 



pi + pi 4 ? 



4q 2 -p't+Po Pt 



for p 2 < p 2 , 



for p\ > pi, 



(3) 
(4) 



with the infrared cutoff parameter p^ = 0.5 GeV 2 . 

Part of the t (i.e. p 2 ) dependence of the amplitude is given by the bottom part of the graphs 
of fig. H], i.e. by the spatial distribution of the colour charge inside the proton. We parametrise 
this distribution by means of the two-gluon form factor F^ G (t) (with F^ G (0) = 1). This form 
factor cannot be calculated within pQCD and should be measured experimentally. The simplest 
hypothesis is that the function Fff*(t) is close to the electromagnetic form factor and can thus 
be written, for instance, as [F^(t)] 2 = e Bt , with B w 5 GeV~ 2 0; equivalently the dipole 
approximation can be used. This hypothesis is implemented in the program. Unfortunately, 
within the LLA there is no shrinkage of the diffractive cone. Therefore the slope B does not 

2 However estimates based upon QCD sum rules p6l give somewhat smaller values of B. 
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depend on energy. This is a shortcoming of the model, and to imitate the effect of the diffractive 
cone shrinkage one should choose different values of B for different energies. 

The remaining part of the t-dependence is contained in the factor {f /[q 2 (2q 2 — pf) In (8q 2 /pl)]} 2 , 
which comes from the upper fermion loop in fig. as well as in the x and q 2 dependence of the 
gluon distribution. 

As mentioned above, a non-relativistic meson wavefunction is assumed. It is possible however 
to allow for relativistic effects by shifting the overall normalisation by a factor rj 2 . For p° mesons 



7] ~ 1.8 |2jJ, while it is expected to be close to unity for J/ip mesons. 



There are a few restrictions on the applicability of the model. In the first place, to be sure 
that |t m m| < Aqcd ( or l^minl < with Rn the nucleon radius), only the region x < 0.1 

should be considered. 

Furthermore, strictly speaking, one should use the pQCD result only at sufficiently large 
values of q 2 , say q 2 > 2 GeV 2 . Nevertheless, the program can be used also at smaller values of q 2 , 
e.g. for p° photoproduction, for which q 2 ~ 0.15 GeV 2 . This is the reason for modifying formula 
(6) of ref. |n|] into (0) by including an infrared cutoff. There are no singularities for q 2 > Aq CD ; 
however the applicability of pQCD for these values of q 2 is questionable. Parametrisations 
can be found for xg(x,q 2 ), with which the experimental results on p° photoproduction are 
reproduced [ ]IDL [ITJ; in this region however it is not clear whether the function xg(x, q 2 ) can still 
be interpreted as the gluon momentum density in the proton. 

We have so far discussed only the case of transversely polarised photons. The exchange of a 
pQCD (two-gluon) pomeron conserves s-channel helicity and, if the vector meson wavefunction 
has the non-relativistic form, the longitudinal cross section is given by a L = (Q 2 /my)cr T . In 
the program both the a L and a T terms are taken into account with the corresponding angular 
distributions for the vector meson decay products. As an example, in the case of two-body 
decays, the cross section is written as: 



d 2 a T 1 

— a T kx (cos ip), (5) 



d cos ddif) An 

-a L k L {costi,?fj), (6) 



d 2 a L 1 



d cos §dip An 

where $ and ip are defined in table |2] and %p = ip + 5, with 5 the angle that the photon 
polarisation vector forms with the normal to lepton scattering plane; this angle is taken to be 
zero with probability (2 — y) 2 /Ny 2 and 90° with probability 1/N, where iV = 1 + (2 — y) 2 /y 2 . 

The angle ip is the difference of the azimuthal angles ip and $, which are defined in two 
different frames. These frames however differ by a boost which is very nearly (to an accuracy 
better than (Q 2 + p 2 )/s) perpendicular to the quantisation axes of both frames, thus hardly 
affecting the azimuthal distributions. 

For decays to spin-1/2 particles (e.g. J/tp — > e + e~), one has /cr(cosf?, ip) = 3/2(1 — 
sin 2 °& cos 2 and k^cos'd) = 3/2 sin 2 For decays to spin-0 particles (e.g. p° — > 7r + 7r~), 

3 There was a mistake in sect. 3 of ref. pl| , where the helicity flip amplitude Ax,o wa s indicated to be non-zero. 
One of us (M.R.) would like to thank D. Kriicker for pointing this out. 



6 



the angular distributions become /^(cost?, if)) = 3 sin 2 $ cos 2 if> and ^(cosi?) = 3(1 — sin 2 ??), 
respectively. 

In the case of the three-body decay V — ► 7r + 7r~7r°, the angular distributions are given by 
fcr(cosi?, "0) oc sin 2 cos 2 -0 and ^(cosi?) oc (1 — sin 2 ??), where now $ and -0 are defined as in 
the two-body decay case but for the normal to the decay plane rather than for the 7r + direction. 

For all other decay modes, the function k is taken as constant. 

The 7p cross section is related to the ep (or pp) cross section as follows: 



dMe d P y^ PV) = TTaT{lP VP) + VL(jL{lP Vp) > (7) 

where Y T and T L are the fluxes of transversely and longitudinally polarised virtual photons, 
respectively. Their explicit expression is: 



2ttQ 2 



l + (l-y) 2 2(l-y)Q 



2 

min 



y Q 2 

ft 2(1 -y) 



(9) 



2ttQ 2 y ' 

where Q 2 i n = M 2 -^ is the minimum value of Q 2 kinematically allowed. 

3 The program 

DIPSI generates events according to reaction (|l|) in which V can be, at present: 

• p° — > 7T + 7r~, 

• UJ — > 7T + 7T _ 7r or 7T + 7T~, 

• </> -> ir+iT", KlKg Or TT+TT-TT , 

• p(1450) -> 7T + 7T~, 7T + 7T _ p Or 7T 7T p , 

• p(1700) -> 7T + 7T-, TT + TT-p OT 7r°7r°p , 

• J/^ — > e + e~, p + p~ or 7r + 7r~7r°, 

• ^'(3600) -> e+e", p+p", ir+ir-J/if) or ttVJ/^, 

• T — ► e + e~, p + p~ or 7r + 7r~7r°. 

The program is steered by control cards, in which the user should fix, for each run, some 
parameters and choose among a few options: 
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• lepton beam type (electron, muon); 

• beam momenta (for fixed target kinematics the proton beam momentum has to be set to 
0; the lepton beam momentum must always be negative); 

• generation limits for Q 2 , y and pf and type of generation for each of these variables (as 
described in section [TT]) ; 

• type of meson produced, decay mode and mass spectrum; 

• value of a s , gluon density and proton form factor; 

• amount of output data. 

Before discussing the structure of the program and its features, we describe the method of 
weighted events on which the generator is based. 



3.1 The method of weighted events 



Each event generated with DIPSI is assigned a weight equal to its cross section. This makes it 
possible to decouple the generation of the independent variables from their expected distribution. 

As an example, let us consider the p\ distribution. The model predicts an approximately 
exponential distribution of the type Aexp (—bp 2 ); suppose that for a given vector meson b = 
10 GeV -2 . If events are generated according to such a distribution it may take considerable 
computer time to substantially populate the high p 2 t region. In order to study the high pj tail, 
one can instead choose to generate with, say, 6 = 3 GeV -2 . To obtain a physically meaningful 
p 2 distribution one has then to weight each event with its cross section and with a factor, that 
we will name "phase space factor" , which takes into account the way in which the variable has 
been generated. 

Assume, for the sake of simplicity, that the model had only one independent variable, £. 
Suppose that, given a random number Ri between zero and unity, sampled from a uniform 
distribution, one obtains a value £j for £ using some algorithm £ = which is able to 

populate adequately the region of £ under study. Based on the model, the value of da/d^i, the 
differential cross section da/d^ for £ = £ i; can then be computed. The weight of the event is 
given by 

da d^ 

= 565ft' (10) 



where the phase space factor d^/dRi is the value of d^/dR computed for R = R^. The phase 
space factor would be unity if £ were sampled from a uniform distribution in the range 0-1. 

The cross section for the process, integrated over the generation range of £, is thus 



a 



N to t i= i 



E (ii) 
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where N tot is the total number of generated events. The statistical uncertainty on a is given by 



A 



a 



N to t\ 



Ntot 
i=l 



(12) 



Likewise, the differential cross section da/d^ for £ = & is obtained as 

1 



tot 



(13) 



where the sum is extended to the events in a narrow bin of width A£ centred around £j. 



3.2 The structure of the program 

The code consists of three main routines: DIPSINI, which controls the initialisation phase; 
DIPSIGEN, for the actual event generation; DIPSOUT, which terminates the program. 



3.2.1 Initialisation 

In DIPSINI the relevant variables are initialised, the random number seed and the external 
input cards are read in and the histograms and n-tuples are booked. The input quantities read 
from the control cards are described in appendix 0. 



3.2.2 Event generation 

For each event the following independent kinematic variables are generated (subroutine DIPSI- 
GEN): y, Q 2 , the mass of the vector meson my, p\ and the azimuthal angles of the scattered 
lepton with respect to the incoming beam direction and of the vector meson with respect to the 
virtual photon direction. The four-momenta of the scattered lepton, of the vector meson and 
of the scattered proton are then fixed by energy-momentum conservation. 

Various options, detailed in appendix 0, are available for generating y (flat, l/y), Q 2 (flat, 
l/Q 2 , 1/Q 4 ) and my (flat, Breit-Wigner, Soding [pR ITU], relativistic p-wave Breit- Wigner) . The 



variable p t is sampled from an exponential distribution of which the user can choose the slope. 

Once the kinematics of the event has been generated, the cross section is evaluated as dis- 
cussed in sect. ^| by subroutine JPRYSKIN, called by DIPSIGEN. Immediately after calling 
JPRYSKIN, DIPSIGEN evaluates the weight by multiplying the cross section by the phase 
space factors. In addition to the weight for the process ep — > epV, that for the reaction 
7*p — > Vp is also computed. 

Finally subroutines JDK (for two-body decays) and TREDK (for three-body decays), also 
called by DIPSIGEN, let the vector meson decay. At this point the weight is modified so as 
to take into account the dependence of the cross section on the decay angular variables; these 
variables are evaluated by the routines DPHELI and HELOMEGA for two-body and three-body 
decays, respectively. 
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3.2.3 Termination 



Subroutine DIPSOUT computes the integrated ep and 7*p cross sections according to expres- 
sion (|TTD and closes histogram and n-tuple files. 



3.3 Implementation and usage 



DIPSI is written in standard FORTRAN 77 and uses the PATCHY offline editor J29| to allow 
easier modifications. All routines are contained in a PAMfile; a short "cradle" program contains 
the list of the routines that the user wants to include and the modifications to the code, if 
any, that need to be applied. By running PATCHY on the PAMfile a FORTRAN program is 
produced. 



The CERN libraries PACKLIB and KERNLIB [ 30 1 are necessary; extensive use is made of 



the format free reading package FFREAD |3T[ and of the histogramming package HBOOK [32| 



(which is a part of PACKLIB). The parton distribution package PDFLIB |33fl may be needed 
(cf. appendix p.l| ). 



3.3.1 Input 

All necessary input parameters are read in from the user control cards (logical unit 8) described 
in appendices |B| and A second input file (logical unit 9) is used in order to keep track of the 
random number seeds. These are the only input files. 

Among the input quantities are the gluon distribution, the strong coupling constant and 
the proton form factor. Concerning these parameters a remark is in order. The cross section 
factorises in a 2 ,, in the gluon density squared and in the proton form factor squared. It is 
thus not necessary to generate new samples of events if one exists already and if any of these 
quantities is modified: it is sufficient to reweight the events. As an example, if events have 
been generated with a gluon density xg{x) = 3(1 — x) 5 and one wants to change to a different 
distribution xg'(x,q 2 ), the weight of each event should be multiplied by the factor [(xg')/(xg)] 2 . 
If desired, however, a user provided parametrisation can be included in subroutine USRGLS; in 
that case the switch USRGLU in the control cards should be set to 1. The default version of 
the subroutine contains, as an example, a call to the MRSA parametrisation from the library 
PDFLIB |[j3||: the user can replace it with any other. 



3.3.2 Output 

DIPSI produces the following output files: 

1. DIPSI. TST (logical unit 18), a log-file listing the values of the input parameters selected 
via the control cards, the integrated ep (or /j,p) and the integrated j*p cross sections; 

2. DIPSI. OUT (logical unit 6), containing the same log-file along with the HBOOK his- 
tograms in ASCII format; 
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Table 3: Cards used to produce the sample of events with elastic production of p° mesons in 
the range 4 < Q 2 < 100 GeV 2 at HERA, cf. figs. |3] and ^. See appendix C for an explanation 
of the meaning of the cards. 



3. DIPSI.HROUT (logical unit 40), with the same histograms readable with PAW [M 

4. DIPSI.NTP (logical unit 41), containing the n-tuple; 

5. DIPSI.ERR (logical unit 19), containing (self-explanatory) error messages, if any; 

6. DIPSI.UNO (logical unit 17), containing all kinematic variables and momenta for a single 
event. The event number is chosen via the card JEVE. 



Figures || to |] show some typical plots. Figures § and [| show the Q 2 , y, W, t and decay 
angles distributions for a sample of p° mesons generated in the HERA kinematic regime, with 
4 < Q 2 < 100 GeV 2 and 0.01 < y < 0.99. Figures || and | show the same plots for a J/ip 
sample produced in the photoproduction regime at HERA. For this sample Q 2 < 4 GeV 2 and 
0.01 < y < 0.99. Tables || and § list the cards used to produce these samples (see appendix 
for a description of the cards). 



4 Summary 

We presented a Monte Carlo generator for elastic vector meson production via reaction ([!]), 



based on a leading logarithm QCD model calculation |2T] 



The model parameters are the strong coupling constant, the two-gluon proton form factor 
and the gluon distribution in the proton. The cross section is proportional to the square of the 
gluon momentum distribution in the nucleon. 



11 



FRF A M 


97 ^ 


TDT<"T FP 


n 
u 


PRF A M 


&9P) no 


TM A SJPF 


n 
u 




1 
1 




n 
u 




n 
u 


M A <3M A Y 


n 
u 




i nnnnn 

1UUUUU 


TTSRC1T ,TT 


n 


TFVF 


n 
u 




n 
u 


VTV/fTNT 
I IvlllN 


1 TT 9 




n 
u 




n go 


TFORFA 


n 
u 


YGEN 


n 


FORFAS 


9 so 


QSQLOW 


l.E-12 


BIPT 


3.00 


QSQUP 


4. 


PTMIN 


0.000 


KEWGEN 





PTMAX 


10.000 


JMESON 





ALPHAS 


0.25 






ETA 


1. 



Table 4: Cards used to produce the sample of events with elastic photoproduction of J/i[) mesons 
at HERA, cf. figs. |5| and || See appendix C for an explanation of the meaning of the cards. 
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Figure 3: Distributions over Q 2 , y, W and t for events with elastic production of p° mesons in 
the range 4 < Q 2 < 100 GeV 2 at HERA. See table | for the list of the cards used. 
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Figure 4: Distributions over cost?, (p, $ and for events with elastic production of p° mesons 
in the range 4 < Q 2 < 100 GeV 2 at HERA. See table || for the list of the cards used. 
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Figure 5: Distributions over Q 2 , y, W and t for events with elastic photoproduction of J/ip 
mesons at HERA. See table £| for the list of the cards used. 
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Figure 6: Distributions over cos??, ip, $ and ^ for events with elastic photoproduction of J /if) 
mesons at HERA. See table £| for the list of the cards used. 
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The program correctly describes the ZEUS data on photoproduction of p° [10, 11] and J /if) 
mesons jyj, and has been used to evaluate the ZEUS acceptance for these and other processes. 
It is interesting that, with suitable input for the gluon distribution, DIPSI reproduces also the 
p° photoproduction data, which are in a q 2 region (q 2 w 0.15 GeV 2 ) below that of applicability 
of perturbative QCD. 

DIPSI is steered by control cards by means of which the user can select, among other features, 
the type and decay mode of the vector meson, the kinematic regime, the range of the independent 
variables and the distribution from which they are sampled. 

Each event is assigned a weight proportional to its cross section. It is thus possible to generate 
the independent variables according to distributions different from those expected by the model, 
thereby enhancing, if necessary, the statistical significance of the sample in regions where the 
cross section is small. 
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A List of the files provided 

The following files are provided: 

DIPSI24.CAR 
DIPSI. CRA 
DIPSI. CARDS 
DIPSI.COM 
DIPSI. UNIX. 



They are respectively the PATCHY pamfile (card version) and cradle, the control cards, and 
the command-files in VMS and in UNIX. The cradle file includes a PATCHY switch to be set 
depending on whether the program is run in a UNIX or a VMS environment. 
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B Selection of model parameters 



B.l The gluon momentum density 

In order to compute the gluon distribution in the proton, two subroutines are provided, GLU- 
ONS and USRGLS. 

In GLUONS a few distributions are implemented. The default one is xg(x) = 3(1 — a;) 5 ; also 
available is the parametrisation obtained by the NMC experiment |3*5| . 

The card ICRXGX (cf. appendix £]) selects the available parametrisations: 

• ICRXGX=0: xg(x) = 3(1 - x) 5 ; 

• ICRXGX=1: NMC lower limit for f = 7 GeV 2 ; 

• ICRXGX=2: NMC central value for f = 7 GeV 2 ; 

• ICRXGX=3: NMC upper limit for q 2 = 7 GeV 2 . 

Furthermore, if ICRXGX=2: 

1. IQ2EVO=l: a q 2 dependent parametrisation |36| of the NMC results is used (fig. |7]); 

2. IQ2EVO=2: same as IQ2EVO=l but xg{x, q 2 ) = xg(x, q^J for q 2 < q 2 mhl = 2 GeV 2 ; 

3. IQ2EVO=3: same as IQ2EVO=2 but xg(x, q 2 ) = x m i n g(x m i n , q 2 ) for x < x min = 3 x 10~ 3 . 

The routine USRGLS is a dummy routine with arguments x, q 2 and xg which returns the 
value of xg(x,q 2 ) for the given x and q 2 . The user can implement the desired parametrisation 
here. The default version of the subroutine contains, as an example, a call to the MRSA 



parametrisation from the library PDFLIB |33 



Finally it is possible to switch from GLUONS to USRGLS by means of the integer card 
USRGLU: USRGLU = 1 allows the use of USRGLS (see appendix 0). 



B.2 The two-gluon form factor of the proton 

The choice of the proton form factor can be done directly via cards (cf. appendix y). The card 
IFORFA allows to choose between the dipole form factor, F N (t) = 1/(1 - g^) 2 , (IFORFA=0) 
and an exponential one (IFORFA=l). In the latter case the slope of the exponential can be set 
via card FORFAS. 



B.3 The strong coupling constant a s 

The value of the strong coupling constant a s is determined by card ALPHAS (see appendix 0). 
A card value in the range between and 1 fixes a s to the card value. ALPHAS values outside 
the range 0-1 let a s evolve with q 2 , assuming a s (q 2 ) = 127r/[25 In (<? 2 /Aq CD )]; in this case, if 
a s comes out to be larger than 0.7, it is automatically set to 0.7. The value of A.qcd is set to 
200 MeV. 
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Figure 7: The gluon distribution in the proton as parametrised in 



B.4 Relativistic effects in the meson wavefunctions 

It is possible to shift the normalisation of the cross section to account for relativistic effects in 
the meson wavefunctions. The cross section can be multiplied by a factor [rj 2 ]; the user can 
choose it by setting the card ETA to the desired value. Zero or negative values are equivalent 
to setting rj = 1. For the p° meson a value 77 = 1.8 is recommended [pTtl. For J/ip mesons rj ~ 1. 



C Input control cards 

An example of input control card file (logical unit 8), with a description of each entry, is given 
below. For more details on the cards controlling the gluon distribution, the form factor of the 
proton, a s and 77, see appendix |TJ. 

C *** Lepton beam momentum in GeV (must be negative...) 
EBEAM -27 . 5 

C *** Proton beam momentum in GeV 
PBEAM 820 . 00 

C *** ntpflag = 0/1: do not fill/fill n-tuple 
NTPFLAG 1 

C *** emc = 0: electron beam, emc = 1: muon beam 
EMC 

C *** Total number of events to be generated 
NUT0 11000 
C *** Event to be dumped in dipsi.uno 
JEVE 4 

C *** y generation range and spectrum: YGEN = — > 1/Y, 
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C *** YGEN = 1 — > flat 

YMIN 4.E-2 

YMAX 9.E-2 

YGEN 1 

C *** Q2 generation range (in GeV2) and spectrum: 

C *** KEWGEN = — > 1/Q2, KEWGEN = 1 — > 1/Q4, KEWGEN = 2 — > flat 

QSQLOW . 10E-9 

QSQUP 5 . 00E+0 

KEWGEN 

C *** JMESON: vector meson type 

C J/psi (also any integer other than the ones below) 

C 1 rho 

C 2 phi 

C 3 Upsilon 

C 4 omega 

C 10 psi' (3686) 

C 11 rho' (1450) 

C 21 rho' (1700) 

JMESON 

C *** JDKLEP : decay mode 

C mu+ mu- (also any integer other than the ones below) 

C 1 e+ e- 

C 2 pi+ pi- 

C 3 K+ K- 

C 4 KOs KOI 

C 5 pi+ pi- piO (with piO — > gamma gamma) 

C 10 pi+ pi- psi — > pi+ pi- mu+ mu- 

C 11 pi+ pi- psi — > pi+ pi- e+ e- 

C 12 pi+ pi- rho — > pi+ pi- pi+ pi- 

C 15 piO piO psi — > piO piO mu+ mu- (piOs do not decay) 

C 16 piO piO psi — > piO piO e+ e- (piOs do not decay) 

C 17 piO piO rho — > piO piO pi+ pi- (piOs do not decay) 

JDKLEP 1 

C *** IMASGE: Vector meson mass distribution 

C Breit Wigner (also any integer other than the 

C ones below) 

C 1 flat 

C 2 Soeding (cf. ref . [28] and [10]) 

C 3 relativistic p-wave BW 

IMASGE 

C *** MASS GENERATION RANGE 

C 

MASMIN 1.00 

MASMAX 2 . 40 

C *** Input gluon distribution (cf. appendix B.l) 

USRGLU 

ICRXGX 
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C *** Q2 evolution: IQ2EVO0 — > no q2 evolution, 



C *** IQ2EV0=1 — > NMC xg(x,q2) 

C *** (cf. appendix B.l) 
IQ2EV0 

C *** Proton form factor: IF0RFA=0 — > dipole, 
C *** IF0RFA=1 — > exponential 

C *** ( c f . appendix B.2) 
IFORFA 



C *** Form factor slope (cf . appendix B.2) 
FORFAS 2 . 50 

C *** pt2 (exponential) generation spectrum and range (in GeV2) 
C *** BIPT = generation slope (in GeV-2) 
PTMIN 0.000 
PTMAX 1000 . 000 
BIPT 5.00 

C *** Alpha strong (cf . appendix B.3) 
ALPHAS . 25 

C *** normalisation shift to account for 

C *** relativistic wavef unction (cf . appendix B.4) 

ETA 1. 

C 

STOP 



D Contents of the n-tuple 

Logical unit 41 is an unformatted file in which we store many of the relevant kinematic variables 
as an n-tuple. In the following we give a brief description of the contents of this n-tuple. For a 
definition of the symbols the reader is referred to tables [I] and §. The n-tuple ID number is 666. 

The components of particle momenta are given in the following order: x, y and z components 
of three-momentum, energy and mass. The z axis coincides with the incoming proton direction 
and is opposite to the direction of the incoming electron. If the proton is at rest the z axis 
is also taken to be opposite to the direction of the incoming electron. Energy, momenta and 
masses are expressed in GeV; cross sections are in nb. Lines from 56 to 70 are relevant only for 
three-body decays. For these decays $ and <p are the polar and azimuthal angles of the normal 
to the decay plane in the helicity frame. The term 'unstable meson' refers to tt°, p° and J/ip. 

* NTUPLE ID= 666 ENTRIES= 100000 DIFFRACTIVE JPSI EVENTS 



* Var numb * Name * Lower * Upper * 



* 


1 


* 


Q2 


* 


0. 


. 100249E-09 


* 


0. 


.499924E+01 


* 


Q 2 


* 


2 


* 


Y 


* 


0. 


.100000E-01 


* 


0. 


.499995E+00 


* 


y 


* 


3 


* 


NU 


* 


0. 


.480671E+03 


* 


0. 


, 240333E+05 


* 


V 


* 


4 


* 


PT2CM 


* 


0. 


, 194709E-05 


* 


0. 


.466549E+01 


* 


pi 


* 


5 


* 


WSQ 


* 


0. 


899670E+03 


* 


0. 


.451005E+05 


* 


w 2 
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* 


6 


* 


Z 


* 


0. 


.997422E+00 


* 


0. 


. 100000E+01 


* 


z 


* 


7 


* 


T 


* 


— 


.600013E-03 


* 


0. 


.467034E+01 


* 


-t 


* 


8 


* 


XL 


* 


0. 


.979781E+00 


* 


0. 


.999787E+00 


* 




* 


9 


* 


PT 


* 


0. 


. 139538E-02 


* 


0. 


.215999E+01 


* 


Pt 


* 


10 


* 


XBAR 


* 


0. 


.212857E-03 


* 


0. 


. 168337E-01 


* 


X 


* 


11 


* 


Q2BAR 


* 


0. 


.239770E+01 


* 


0. 


.423783E+01 


* 


f 


* 


12 


* 


HCOSTH 


* 


- 


.999990E+00 


* 


0. 


.999992E+00 


* 


cos-# 


* 


13 


* 


HPHI 


* 


0. 


. 123719E-02 


* 


0, 


.628306E+01 


* 




* 


14 


* 


HPHIC 


* 


0. 


. 184483E-02 


* 


0. 


.628304E+01 


* 


$ 


* 


15 


* 


HPS I 


* 


0. 


. 150895E-02 


* 


0. 


.628255E+01 


* 


i> 


* 


16 


* 


WEIGHT 


* 





.897848E-05 


* 


0. 


.546283E+01 


* 


ep weight 


* 


17 


* 


WTGAMP 


* 





.565991E-04 


* 


0. 


.263028E+02 


* 


7P weight 


* 


18 


* 


EBE1 


* 


0. 


.OOOOOOE+OO 


* 


0. 


. 000000E+00 


* 


incoming 


* 


19 


* 


EBE2 


* 


0. 


.OOOOOOE+00 


* 


0. 


. 000000E+00 


* 


lepton 


* 


20 


* 


EBE3 


* 


— 


.275000E+02 


* 


— 


. 275000E+02 


* 




* 


21 


* 


EBE4 


* 


0. 


.275000E+02 


* 


0. 


. 275000E+02 


* 




* 


22 


* 


EBP1 


* 


0. 


.OOOOOOE+OO 


* 


0. 


. OOOOOOE+OO 


* 


target 


* 


23 


* 


EBP 2 


* 


0. 


.000000E+00 


* 


0. 


. OOOOOOE+OO 


* 


proton 


* 


24 


* 


EBP3 


* 


0. 


.820000E+03 


* 


0. 


. 820000E+03 


* 




* 


25 


* 


EBP4 


* 


0. 


.820001E+03 


* 


0. 


.820001E+03 


* 




* 


26 


* 


ESE1 


* 


— 


.221037E+01 


* 


0. 


.219890E+01 


* 


scattered 


* 


27 


* 


ESE2 


* 


- 


.220185E+01 


* 


0. 


.219534E+01 


* 


lepton 


* 


28 


* 


ESE3 


* 


— 


.272250E+02 


* 


- _ 


. 137196E+02 


* 




* 


29 


* 


ESE4 


* 


0. 


. 137501E+02 


* 


0. 


.272669E+02 


* 




* 


30 


* 


ESE5 


* 


0. 


.511000E-03 


* 


0, 


.511000E-03 


* 




* 


31 


* 


ESP1 


* 


— 


. 179924E+01 


* 


0. 


. 161310E+01 


* 


recoiling 


* 


32 


* 


ESP2 


* 


— 


. 197851E+01 


* 


0, 


. 187068E+01 


* 


proton 


* 


33 


* 


ESP3 


* 


0. 


.803420E+03 


* 


0, 


.819825E+03 


* 




* 


34 


* 


ESP4 


* 


0. 


.803421E+03 


* 


0. 


.819826E+03 


* 




* 


35 


* 


ESP5 


* 


0. 


.938272E+00 


* 


0. 


.938272E+00 


* 




* 


36 


* 


GAM1 


* 


— 


.219890E+01 


* 


0. 


.221037E+01 


* 


virtual 


* 


37 


* 


GAM2 


* 


- 


.219534E+01 


* 


0. 


.220185E+01 


* 


photon 


* 


38 


* 


GAM3 


* 


- 


. 137804E+02 


* 


- _ 


. 275000E+00 


* 




* 


39 


* 


GAM4 


* 


0. 


.233082E+00 


* 


0. 


. 137499E+02 


* 




* 


40 


* 


GAM5 


* 


— 


.223590E+01 


* 


0. 


.428417E-02 


* 




* 


41 


* 


VEC1 


* 


— 


.306449E+01 


* 


0. 


.294946E+01 


* 


vector 


* 


42 


* 


VEC2 


* 


— 


.298695E+01 


* 


0. 


.297885E+01 


* 


meson 


* 


43 


* 


VEC3 


* 


- 


. 135730E+02 


* 


0. 


. 162611E+02 


* 




* 


44 


* 


VEC4 


* 


0. 


.309709E+01 


* 


0. 


. 168128E+02 


* 




* 


45 


* 


VEC5 


* 


0. 


.309690E+01 


* 


0. 


.309690E+01 


* 




* 


46 


* 


MUP1 


* 


— 


.322182E+01 


* 


0. 


.326902E+01 


* 


first 


* 


47 


* 


MUP2 


* 


- 


.339305E+01 


* 


0. 


.353246E+01 


* 


decay 


* 


48 


* 


MUP3 


* 




. 136807E+02 


* 


0. 


.127291E+02 


* 


particle 


* 


49 


* 


MUP4 


* 


0. 


. 182204E+00 


* 


0. 


, 136921E+02 


* 




* 


50 


* 


MUP5 


* 


0. 


.511000E-03 


* 


0. 


.511000E-03 


* 




* 


51 


* 


MUM1 


* 




.303850E+01 


* 


0, 


.310788E+01 


* 


second 


* 


52 


* 


MUM2 


* 




.338481E+01 


* 


0, 


.358694E+01 


* 


decay 
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* 


53 


* 


MUM3 


* 


- , 


. 135949E+02 


* 


0. 


. 127755E+02 


* 


particle 


* 


54 


* 


MUM4 


* 


0. 


. 188318E+00 


* 


0. 


. 136271E+02 


* 




* 


55 


* 


MUM5 


* 


0. 


.511000E-03 


* 


0. 


.511000E-03 


* 




* 


56 


* 


PIZ1 


* 


0. 


OOOOOOE+00 


* 


0, 


. OOOOOOE+00 


* 


third decay 


* 


57 


* 


PIZ2 


* 


0. 


OOOOOOE+00 


* 


0, 


. OOOOOOE+00 


* 


particle 


* 


58 


* 


PIZ3 


* 


0. 


OOOOOOE+OO 


* 


0, 


. OOOOOOE+00 


* 


( "unstable" 


* 


59 


* 


PIZ4 


* 


0. 


.OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


meson: 7r°, p°, 


* 


60 


* 


PIZ5 


* 


0. 


OOOOOOE+OO 


* 


0, 


. OOOOOOE+00 


* 




* 


61 


* 


PHI 


* 


0. 


.OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


first 


* 


62 


* 


PH2 


* 


0. 


.OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


particle from 


* 


63 


* 


PH3 


* 


0. 


OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


"unstable" 


* 


64 


* 


PH4 


* 


0. 


.OOOOOOE+00 


* 


0. 


.OOOOOOE+00 


* 


meson 


* 


65 


* 


PH5 


* 


0. 


.OOOOOOE+00 


* 


0. 


.OOOOOOE+00 


* 




* 


66 


* 


PH01 


* 


0. 


OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


second 


* 


67 


* 


PH02 


* 


0. 


OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


particle from 


* 


68 


* 


PH03 


* 


0. 


.OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


"unstable" 


* 


69 


* 


PH04 


* 


0. 


OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 


meson 


* 


70 


* 


PH05 


* 


0. 


OOOOOOE+00 


* 


0. 


. OOOOOOE+00 


* 





E List of routines and functions 



BREIT 


non-relativistic Breit-Wigner mass distribution 


DCR0SS 


vector product between two vectors 


DEPHI 


azimuthal angle determination 


DIPSI 


main program 


DIPSIGEN 


generation of kinematic variables 


DIPSINI 


initialisation 


DIPSOUT 


cross section evaluation and output 


DPHELI 


helicity angles determination 


GLUONS 


evaluation of gluon distribution 


JDK 


two-body decay 


JPRYSKIN 


cross section determination 


HELOMEGA 


helicity angles for three-body decays 


LORENZ 


three dimensional Lorentz boost 


LUBEND 


closing of unformatted file (lun 40) 


LUBOOK 


histogram and n-tuple booking 


LUFIL 


histogram and n-tuple filling 


NTPINIT 


definition of n-tuple variables 


RUOTA 


reference system rotation 


SODING 


Soeding shape distribution 


SUB 5 


subtraction of two vectors 


TIMTEM 


time usage evaluation 


TREDK 


three-body decay (including decay of most unstable meson) 
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ULANGL 
USRGLS 
WIGNER 



angle determination 

user routine to evaluate the gluon distribution 
relativistic Breit-Wigner mass distribution 
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